Molecular ordering of glycine on Cu(lOO): the p(2 x 4) superstructure 
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Glycine molecules deposited on Cu(lOO) surlace give rise to an anisotropic free-electron-like (FEL) 
electronic dispersion in its p(2x4) superstructure, as reported in recent experiments [Phys. Rev. 
Lett. 99, 216102 (2007); J. Am. Chem. Soc. 129, 740 (2007)]. Using density functional theory and 
exhaustively calculating sixteen possible structures, we have determined the molecular arrangement 
that can give the experimentally observed FEL behavior. Eight configurations, among the sixteen, 
were not investigated before in the literature and one of them (denoted Str-3) is able to provide the 
FEL behavior in excellent agreement with the experiments. In addition, the particular configuration 
Str-3 satisfies other criteria of the observed p(2x4) superstructure, e.g. chirality and cleavable 
orientation. 

PACS numbers: 81.07.Nb 68.43.Hn, 73.20.At 



I. INTRODUCTION 

It was demonstrated that nearly free electron surface 
states in low-index surfaceJ^' of noble metal, e.g. Cu, 
Ag, and Au, can be tuned by surface adsorbates^El^ step 
edgetP^, artificial nanostructure^, confinement along 
certain dimension^, and among the others. Molecules, 
due to its unique capability of building various nanos- 
tructures, e.g. nano-gratings-- and nano-mest^, were 
employed to manipulate the surface electronic states of 
metals. Recently, free-electron-like (FEL) behavior was 
discovered on the molecular nanostructurcs supported by 
metal surfaces, rather than on the metal surfaces them- 
selves, which has attracted great interest in both physics 
and chemistry communitie^iitt^Sl xhe FEL states appear 
to show different electron effective masses depending on 
the material detail, demonstrating that surface electronic 
structure can be engineered t o produce desired carrier 
mobility in molecular system^SHSIIlIini, As an exam- 
ple, a FEL behavior was observed in perylene-3,4,9,10- 
tetracarboxylic-3,4,9,10-dianhydride (PTCDA) molecu- 
lar islands on Ag(lll)P^, which was explained as a con- 
fined 2D interface state hybridized by the lowest unoccu- 
pied molecular state of PTCDA and the surface state of 
Ag underneatlPH2ni_ 

More recently, an anisotropic FEL behavior in a p(2x4) 
superstructure of glycine molecules (in the form o f gly- 
cinate) on isotropic Cu(lOO) surface was reportecP^^^. 
This suggests a potential ability to modulate electronic 
structure and especially the electron effective mass of 
molecular monolayers or molecule-metal interfaces by us- 
ing organic molecules. A key question concerning this 
interesting molecular surface is why electrons move with 
different speeds in two different directions. In other 
words, why an isotropic Cu(lOO) surface produces an 
anisotropic FEL behavior after glycine molecules are de- 
posited. Clearly, the geometric structure of the molecu- 
lar system plays an essential role in generating such an 
anisotropic FEL behavior. However, so far the issue of 



geometric structure has not been adequately addressed in 
the literature. It is the purpose of this paper to report a 
systematic investigation which reveals the proper molec- 
ular configuration that supports the anisotropic FEL dis- 
persion. 

About a decade ago, scanning tunneling microscopy 
(STM) was employed to investigate glycine/Cu(100) by 
Yang & Zhao's group^^IHlI]^ Two superstructures were 
found and their atomic structures were proposed on the 
basis of the STM images and an understanding of hy- 
drogen bonds. One is a c(2x4) homochiral superstruc- 
ture (denoted Str-g), as shown in Fig[l|a). The edges 
of its molecular islands on the surface are along the 
< 310 > direction, as schematically indicated by the 
red solid line. Another superstructure is a p(2x4) hete- 
rochiral superstructure (denoted Str-d, see Fig. [l|^b)), 
in which the edges of its molecular islands are along 
the < 100 > direction. Spots in the associated low 
energy electron diffraction (LEED) study indicate that 
these two phases co-exist-^^'^. Profiting from well devel- 
oped f abric ation techniques, the recently observed STM 
\ma.geff^^^^ clearly show that the preferred edge orien- 
tation, i.e. cleavable orientation, of p(2x4) molecular 
islands is also along < 310 >. (Here we define the 
"< 310 > rule" as that the cleavable orientation of molec- 
ular islands is along < 310 >.) The STM result indicates 
that the p(2x4) structure proposed by Yang & Zhao's 
group, in which the molecular islands cannot be cleaved 
in the < 310 > direction, appears to be inconsistent with 
the recent observatiorP^. 

Photoemission diffraction (PhD) technique provides 
vertical atomic information rather than local density of 
states as provided by STM. In a PhD experimenlP^, a 
joint European group proposed another p(2x4) hete- 
rochiral superstructure, i.e. the (2x4)-pg (denoted Str-h) 
shown in Fig. [ijc). It was found that this structure, as 
a predominant structure, may co-exist with the c(2x4). 
Reference [24l also indicates that the (2x4)-pg structure 
satisfies the < 310 > rule, which thus appears to be con- 




FIG. 1. (Color online) Four molecular configurations, (a) to (c) are the configurations previously proposed and investigated, 
i.e. c(2x4)(Str-g in this work, (a)), p(2x4) (Str-d, (b)), and p(2x4)-pg (Str-h, (c)). (d) is a new configuration proposed in this 
work, in which the in-plane molecular orientation is roughly perpendicular to that in (a) to (c). The red solid line in (a) denotes 
[130], a direction included in < 310 > direction. 



sistent with the recent STM observatiorf^. 

Attempting to clarify the situation, Mae and 
Morikawa^ performed a density functional theory (DFT) 
calculation on the t hree s tructures proposed by the two 
experimental groupd^ilSil, Jt was found that the total 
energy of the proposed c(2x4) (Str-g) and p(2x4) (Str- 
d) structures are energetically very close to each other, 
while that of (2x4)-pg (Str-h) is significantly higher. 
Structure (2x4)-pg (Str-h) was thus deemed to be un- 
likely. In other words, the DFT calculation reported in 
Ref i25, poses a challenge to the community that none of 
the propose d p(2x4) st ructures could satisfy all previous 
experimentaP^ESOlHll! and theoreticaP^ results. 

It was recently reported that the electronic structure 
of p(2 X 4) shows an anisotropic FEL behavior as mea- 
sured by scanning tunneling spectroscopy (STS), which 
provides a new avenue for solving the above discussed 
discrepancy. Particularly, the measured electronic struc- 
ture, i.e. the electron effective mass (m*) and energy 
levels at certain fc-points, can be adopted as criteria to 
assess all the proposed molecular configurations, finding 
or verifying the configuration of the observed p(2 x 4) 
superstructure. 

Based on this idea, we therefore relaxed the atomic 
structure and systematically calculated the total energy 
of 16 configurations including eight newly proposed ones 
(see Fig. [2]). For those fully relaxed configurations that 
show particularly strong stability, we calculated their 
electronic structures. The calculated band structure of 
configuration Str-3, proposed in this work (atomic struc- 
ture shown in Fig. [ijd)), is impressively consistent with 
the STS results in both surface directions. It is re- 
markable that the geometry of Str-3 is highly consistent 
with both the < 310 > rule and the observed chiral- 
ity. Configuration Str-3 is also preferred in the theory- 
experiment comparison of STM images. All these consis- 
tencies with experimental observations that other config- 
urations cannot offer, strongly suggest that Str-3 is most 
likely the configuration describing the experimentally ob- 
served hetero-chiral p(2x4) superstructure in the recent 
STM and STS measurements. 
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FIG. 2. (Color online) According to O-O orientation, 0-N 
direction, the relative position of equivalent atoms, and the 
chirality, there are 2* — 16 configurations that should be con- 
sidered, including eight previously considered configurations 
(Str-a to Str-h) and eight newly proposed configurations (Str- 
1 to Str-8). 



The rest of the paper is organized as follows. In Sec- 
tion II, we present computational details including the 
molecular configurations being considered and investi- 
gated. Section III presents results and discussion. The 
last section provides further discussion and a summary. 



II. COMPUTATIONAL DETAILS 

As discussed in Section [ij the main issue is concern- 
ing the p(2 X 4) superstructure. We thus only focus on 
this structure in the rest of the paper. In particular, the 
atomic structure of c(2 x 4) can be easily derived once the 
structure of p(2 x 4) is determined. It was reported that 
the carboxylic acid group (COOH) in glycine molecules is 
dissociated on Cu(lOO), giving glycinate molecules that 
will be discussed in this paper. Notice that the line joint- 
ing the two O atoms in a glycine (the dashed line in Fig. 
[ija)) is roughly parallel to the [110] direction (the shorter 
edge of a (2x4) lattice, solid rectangle. Fig. Ilia)) in all 
previously proposed and studied configurations. Never- 



3 



TABLE I. Calculated total energy and summarized information (Chirality, cleavable orientation, and the type of lattice) of all 
possible configurations. 



Str- 12345678abcdefgh 
Energy (eV) 1.03 0.60 0.32 0.43 0.55 0.83 0.40 0.22 0.80 0.40 0.23 0.00 0.98 0.70 0.04 0.29 

Chirality NYYNNYYNYNN YYN N Y 
Orientation YYYYYNNNYNNNYYY Y 

X(2 x4) pppcppp-pp - pppc p 



theless, there appears no reason why this 0-0 line cannot 
be parallel to the [110] direction (the longer edge of the 
lattice, solid rectangle. Fig. [ija)). Surprisingly, this pos- 
sibility has not been considered before, e.g. the configu- 
ration shown in Fig. [l|d) where the 0-0 line is roughly 
parallel to the longer edge of the lattice. Based on this 
observation that glycine may take another orientation 
forming a p(2x4) overlayer, 16 possible configurations 
were considered, in which eight of them are proposed by 
this work for the first time as shown in Fig. [2j 

All calculations were done using DFT, the general- 
ized gradient approximation for the exchange-correlation 
potentiaPS, the projector augmented wave metho(p3^ 
and a plane-wave basis set as implemented in the Vienna 
ab-initio simulation package^^'. Six layers of Cu atoms 
and a 2x4 supercell, separated by a 10-layer vacuum re- 
gion were employed to model the Cu(lOO) surface. The 
molecule is only put on one side of the slab and a dipole- 
correction was applied. A fc-mesh of 8 x 4 x 1, verified 
byal2x6xl one, was adopted to sample the surface 
Brillouin Zone for both relaxation and total energy cal- 
culations. The energy cut-off for the plane waves is set 
to 400 eV. In structural relaxation, all atoms except for 
the bottom three layers were fully relaxed until the net 
force on every atom is less than 0.02 eV/A. 

III. RESULTS AND DISCUSSION 

A. Molecular configurations 

Table |l] shows the calculated relative total energies of 
all configurations with respect to the lowest energy which 
is the total energy of Str-d. Configurations Str-1 to Str-8 
were not considered before in the literature. Letter "Y" 
("N") in row "Chirality" indicates that the correspond- 
ing configuration is (not) chiral. The next row shows 
whether or not the configuration can be perfectly cleaved 
along < 310 > directions (the < 310 > rule). Here, term 
"perfect cleaving" of a configuration along a certain di- 
rection means that a Cu mono-atomic step along such di- 
rection can be and only be covered by glycine molecules 
in the molecular arrangement of that configuration. Row 
"X(2 X 4)" shows the type of lattice of all configurations, 
in which Str-4 and Str-g have a c(2 x 4) lattice, the lattice 
of Str-8 and Str-c degrades into a (2 x 2) one, and the 



others is p(2 x 4). 

It was found that the relative total energies of Str-8 
and Str-c are well consistent, which implies a good con- 
vergence of relative energies between different molecu- 
lar orientations to less than 0.01 eV in our calculation. 
Since Str-8 and Str-c degrade into a (2 x 2) structure and 
are not observable in all experiments, they are thus ex- 
cluded from further discussion. From a purely energetic 
point of view, Str-d and Str-g are the two suggested con- 
figurations which were also preferred in the literatur^^H^ 
The energy difference between them is 0.04 eV, consistent 
with the previous DFT calculatioiPI. 

Purely energetic aspect aside, there are two additional 
factors which may influence the determination of the ob- 
served configuration, i.e. the chirality and the cleavable 
orientation. Chirality leads to a larger entropy of a chiral 
structure, giving a lower free energy at finite tempera- 
ture. Although Refs. [2Tti23l did not state the preferred 
orientation of the island edges, it is indeed observable 
that most edges of glycine islands shown in the Fig. 1 of 
Ref. [12] are in the < 310 > direction. Whether or not 
the superstructure is cleavable along the < 310 > direc- 
tion, we believe, is of importance for distinguishing the 
configurations. 

According to Table |Ij Str-d has the lowest total en- 
ergy, chirality, but is not cleavable along the < 310 > 
direction. The c(2 x 4) configuration Str-g has the sec- 
ond lowest total energy, no chirality, but can be cleaved 
along the < 310 > direction. These results are con- 
siste nt wi th the experimental observations of Zhao & 
Yan^2il22| On the other hand, Str-h, the PhD experi- 
ment suggested configuration, has chirality and is cleav- 
able along < 310 >. Among those configurations having 
chirality and also cleavable along < 310 >, Str-h has the 
lowest total energy. It also has the third lowest total en- 
ergy among all configurations (Str-8 and Str-c excluded). 
We thus conclude that the PhD suggested configuration 
is quite reasonable. Table |l] shows that the total energies 
of all configurations proposed by this work are higher 
than the previously proposed ones, i.e. Str-d, Str-g and 
Str-h. However, Str-3, a newly proposed configuration 
that has chirality and is cleavable along < 310 >, has a 
total energy very close to Str-h (within 0.03 eV). There- 
fore, it should be considered as a likely candidate for the 
observed p(2 x 4) superstructure. 

Recent experiments suggested that the p(2 x 4) struc- 
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TABLE II. Structural properties obtained from the DFT calculation in comparison to experimental data. Z]v{o)i stands for 
the difference in z direction between N (O) and the first layer Cu underneath; while Z12 (N(0)) means that between first and 
second layer Cu atoms under N(0). dcti_jv(o) shows the bond length of Cu and N(0). ^jv(o) is the angle between the vector 
perpendicular to the surface and the one formed by Cu and N(0). There are two O atoms in each molecule, so that two values 
are listed in the table. 



Str- 


d 


h 


3 


Exp.M 


Zni (A) 


2.07 


2.08 


2.07 


2.04±0.02 


dcu-N (A) 


2.10 


2.08 


2.09 


2.05±0.02 


^12 (N) (A) 


1.89 


1.89 


1.90 


1.80±0.09 


On n 


8.8 


4.6 


9.7 


5±4 


Zoi (A) 


2.06 (2.06) 


2.05 (2.09) 


2.03 (2.13) 


2.02±0.02 


dcu-o (A) 


2.08 (2.18) 


2.10 (2.20) 


2.06 (2.23) 


2.05±0.02 


Z12 (O) (A) 


1.81 (1.81) 


1.81 (1.83) 


1.82 (1.83) 


1.79±0.06 


oo n 


7.2 (18.1) 


12.3 (17.9) 


10.9 (16.1) 


9±2 



ture should have chiraUty, can be perfectly cleaved along 
< 310 >, show an isotropic free-electron-likc behavior of 
certain bands^^ In addition, as a basic rule, the total 
energy of the structure or the kinetic energy barrier to 
reach it, should be as low as possible. Based on our total 
energy calculation, the analysis of chirality and cleavable 
orientation, we selected three configurations Str-d, Str- 
h, and Str-3, for band structure calculations as described 
in section C below. Reasons are the following: (1) Str- 
d has chirality and holds the lowest total energy among 
all configurations; (2) Str-h and Str-3 satisfy both rules 
of chirality and cleavable orientation, and have relatively 
low total energies. 

B. Atomic structures 

Structural properties that the associated experimen- 
tal data are available for configurations Str-d, Str-h, and 
Str-3, are summarized in Tab. |lT] in which the val- 
ues for Str-d are the same as the theoretical result in 
Ref.f^. Theoretical values of each configuration are con- 
sistent with the corresponding values measured by photo- 
diffractiorPS, hence it is difficult to assign the preferred 
configuration according to such experiment. The table 
shows that the differences in distances and angles be- 
tween these three configurations at the single molecular 
level are rather small, implying that the energetic differ- 
ence is likely not primarily due to the tiny structural. 

The intermolecular interaction, i.e. hydrogen bond- 
ing, is therefore involved in discussion. Each configura- 
tion contains four categories of hydrogen bonds including 
three intermolecular bonds and one intramolecular bond, 
as marked in Fig. [ij^b) . The shortest bond length in Str- 
d is 2.00 A which almost approaches the shorter limit of 
a typical hydrogen bond, indicating a stronger strength 
of this bond. Although the difference between the short- 
est and longest bonds in Str-d is somewhat large, i.e. 
0.68 A , all four bonds can be reasonably considered as ac- 



TABLE III. Bond lengths of the four hydrogen bonds, as 
marked in Fig. [l] of configurations Str-d, Str-h, and Str-3. 



Hydrogen bonds: 


1 


2 


3 


4 


Str-d (A) 


2.61 


2.00 


2.68 


2.36 


Str-h (A) 


2.35 


1.91 


1.70 


3.64 


Str-3 (A) 


2.24 


2.36 


2.12 


2.58 



tual hydrogen bonds. In difference from Str-d, the bond 
lengths in Str-h are distributed in a fairly wide range, 
from 1.70 A to 3.64 A. It therefore leads to a higher 
total energy of Str-h compared with Str-d due to the ab- 
sence of a bond, although the interaction strength of the 
1.70 A and 1.91 A bonds are rather strong. In terms 
of Str-3, its bond lengths are very close to each other, in 
the range from 2.12 A to 2.58 A. The quadrangle formed 
by these four bonds is rather similar to a square, indicat- 
ing a good balance of hydrogen bonding strength in both 
directions. It is still inconclusive which configuration is 
the one that was observed in the experiments showing 
an anisotropic FEL behavior, after atomistic structures 
of those configurations were considered. It calls for a 
further investigation to distinguish these configurations. 



C. Band Structures 

We have calculated the band structures of configura- 
tions Str-3, Str-d, and Str-h along the two directions from 
point G (r) to X and Y, respectively, as shown in Fig. 
H The reciprocal lattice vectors of the p(2 x 4) supercell 
is shown in the inset of Fig. |3] (left panel) . The A;- vector 
corresponds to the [110] direction in real space (in 
blue for all panels) and represents the [110] direction 
in real space (in red for all panels). Both colors are con- 
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FIG. 3. (Color online) Calculated band structures along [110] (blue down-half-fiUed circle) and [110] (red up-half-filled circle) 
of Str-d (a), Str-h (b), and Str-3 (c). Red and blue shaded areas present the experimentally measured curves with a broadening 
of 0.04 eV for the clarity of plotting. Both colors are consistent throughout the paper, i.e. blue for [110] and red for [110]. The 
first Brillouin zone together with three k points (G, X, Y) were shown in (a) inset. 



sistent with experimental plots, i.e. red for the [110] and 
blue for the [110] direction. In order to compare with 
the experimental measuremenlfi^lj the band-structure fig- 
ures only show the disper sion relation s fr om F to the 
points away from F by 0. 2 and 0.4 cf, respectively 
(roughly 0.12 in length). Red and blue shaded ar- 
eas in Fig. [3] present the experimentally measured curves 
with a broadening of 0.04 eV. In the experimenlP^, a 
FEL state in the [110] direction (in red, GY in fc-space), 
with an electron effective mass m^*^ 0.6, was observed 
around 0.11 eV above the Epermi at F. Another FEL 
state, starting from the same energy at the F, with an 
rrie*^ 0.06, was detectable only in the [110] direction 
in fc-space, in blue). 

The calculated band-structure of every configuration 
shows two perfect parabolic bands starting from roughly 
-0.05 eV at the F, which are observable in both directions 
(both in red and blue). These two bands are assigned as 
surface states of the bare Cu surface^-'^ at the bottom of 
the slab. It differs within O.OleV for the calculated val- 
ues among the three molecular configurations, indicating 
again that the energy resolution of our calculation is bet- 
ter than 0.01 eV. 

According to the band-structure of Str-d (left panel 
of Fig. |3|, except the state located at -0.05 eV, only 
one band in red shows a FEL dispersion relation with 
nie*^ 0.6, which resides 0.03 eV above the Epermi- An- 
other band, with a nie*^ 0.06, has almost the same eigen- 
value at the F and reaches 0.35 eV at roughly 0.075 A~^, 
although there is a band-crossing induced gap opening 
from 0.05 eV to 0.13 eV (actual gap from 0.09 eV to 
0.11 eV). Even if the gap opening is not a problem for 
STM-detection, the energy of these two bands at the F 
point, in between 0.03 eV and 0.04 eV, still cannot offer 



a comparable result with the experimental observation. 
We therefore can rule out configuration Str-d as that ob- 
served experimentally. 

The result of Str-h, as shown in Fig. [3] (middle panel) , 
is even further away from the experimental result than 
that of Str-d. Two bands hybridized by the Cu surface 
states and the molecular states of glycine are located 
at 0.06 eV above the Epermi- Apparently, it shows an 
anisotropic FEL behavior with a very small gap opening 
in the direction shown in blue. However, the effective 
masses of these two bands, i.e. smaller (larger) mass 
for the red (blue) band, are opposite to the experiment 
(smaller mass for the blue). We therefore have to exclude 
Str-h as a likely candidate for the p(2 x 4) superstructure 
measured by STS. 

Finally, the right panel of Fig. |3] shows the calculated 
band structure of Str-3. It was found that two bands sit- 
ting around 0.13 eV at the F point are capable of show- 
ing the anisotropic FEL behavior. Both bands show an 
impressive consistency with the experimental measure- 
ments, though slight gap openings appear in the blue 
curves, i.e. at around 0.02 (wave number) / 0.17 eV 
(energy) for the blue band (smaller nif,*). In particular, 
the calculated eigenvalues at different fc points of both 
bands match the measured data fairly well in both di- 
rections, notwithstanding the theoretical band bottoms 
are 0.02 eV higher than the experimental value. The 
theory-experiment agreement is quite impressive because 
the small difference, i.e. 0.02 eV, almost reaches the en- 
ergy resolution of STS measurements and DFT calcula- 
tions. 

The electron effective mass of the two bands in the 
two directions were fitted for configuration Str-3. In the 
red direction ([TlO]), the fitted m* is 0.54, which is 0.07 
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Exp. Str-3 Str-d Str-h 

FIG. 4. (Color online) Experimental STM image of the p(2 x 4) superstructure (Fig. 2(e) in Ref. , 13 j (left) and simulated 
STM images for configurations Str-3 (middle left), Str-d(middle right) and Str-h(right), respectively. The sample bias in the 
simulation is -100 mV, the same as that in the experiment. 



smaller than the measured valu^^ - if all energy points 
in the band are used in the fitting. In the experimen- 
tal pape/^, however, the k value of the first data point 
is around 0.04 A~^. Therefore, we re-fit the theoretical 
data by removing energy points for k in between 0.00 
and 0.04 and the obtained m* is 0.60, which is only 
0.01 smaller than the measured value of 0.6^^. In an- 
other direction, the fitted m* is 0.07, which only differs 
by 0.01 from the experimental value of O.OfP^. Notice 
that the fitted effective mass is somewhat sensitive to the 
points near the T poin t which were not given in the ex- 
perimental paper^i21i3| q^^^. results strongly suggest that 
configuration Str-3 is the one showing the STS measured 
anisotropic FEL behavior. 



D. STM simulation 

The simulated STM images of configurations Str-3, 
Str-d, Str-h were compared with an available experi- 
ment. Figure |4] shows an experimental STM image of 
the p(2 X 4) superstructure, together with the simulated 
STM images for Str-3, Str-d and Str-h, respectively. The 
simulation was performed using the Tersoff-Hamann ap- 
proximation. According to the experimental conditions, 
e.g. -100 mV sample biaJl^, the sample charge density 
was integrated from the states that their energy is in the 
range from -100 meV to the Fermi level (OmeV). An iso- 
surface of charge density at 0.002e/A'^ is plotted in each 
panel. 

Each simulated image contains several bright spots 
that are nearly centered around the upward H atoms con- 
nected with C atoms, as indicated by the atomic struc- 
ture overlaid on the image. The image of Str-h shows a 
hexagonal structure, similar to the corresponding result 
reported in Ref. I25L but which is completely different 
from that of the experiment. In terms of the position 
of the bright spots solely, the image of Str-d appears 
similar to that of Str-3. Both of them more or less re- 
produce the experimental observation and can hardly be 
distinguished. Notice that the experimental image also 
presents a ribbon-like feature, which is offered by Str-3 
but not by Str-d. Although the STM image may be dis- 



torted by the subtle interaction between the tip and the 
upward H atom, the consistency of the ribbon-like fea- 
ture between the experiment and the simulation of Str-3 
suggests again that Str-3 is the most likely candidate for 
the p(2 X 4) superstructure observed in Ref. [131 



IV. FURTHER DISCUSSION AND 
CONCLUSION 

There is an interesting issue that the STM/STS ex- 
periments are consistent with that of Str-3 but not with 
that of Str-d, although Str-d is energetically more favored 
than Str-3. Kinetic effect could most likely be responsible 
for the stabilization of configuration Str-3 against Str-d, 
according to a preliminary DFT calculatiorPSl. As a re- 
sult, when the growth/evaporation rate is low, e.g. as 
that in molecular beam epitaxy, kinetic effect dominates. 
Therefore, Str-3 is kinetically favored, resulting in the 
formation of molecular islands in the Str-3 configuration. 

Notice that the atomic arrangement of N and O of Str- 
3 is the same as that of Str-h. It means Str-3 also satisfies 
the PhD experimental observations^'*. Our analysis pre- 
sented above suggests that the experimental preparation 
conditions in Refs. [12 [T3 and [2^ should be similar, but 
they may substantially differ from that of Refs. l2TI - [23l 
thereby leading to two observed p(2x4) superstructures, 
i.e. the kinetically favored Str-3 and the energetically 
favored Str-d. 

In summary, according to various experimental and 
theoretical results, the long-standing discrepancy regard- 
ing the exact p(2x4) structure of glycine on Cu(lOO) 
has been elucidated. The exact configuration of the 
STS measured p(2x4) superstructure is suggested Str- 
3, a new configuration not considered before, in which 
the 0-0 line in a glycine is oriented roughly parallel 
to the longer edge of the p(2x4) supercell. This con- 
figuration has chirality, satisfies the < 310 > rule, and 
also shows a relative good thermal stability. Its calcu- 
lated band structure shows impressive agreements with 
the measured data, while other configurations cannot of- 
fer such agreements. Given the revealed configurations, 
interesting electronic behavior of this system can be fur- 
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ther investigated. Finally, the Str-d configuration pro- 
posed in Ref. [21] and verified by DFT^S,^ jg expected to 
be observable after an annealing treatment applied to a 
low-rate grown sample where Str-3 is predominant. 
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